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Abstract: Heavy fermion alloys at critical doping typically exhibit non- Fermi-liquid 
behavior at low temperatures, including a logarithmic or power law rise in the ratio 
of specific heat to temperature as the temperature is lowered. Anomalous specific 
heat of this type is also observed in a simple class of gravitational dual models that 
exhibit anisotropic scaling with dynamical critical exponent z > 1. 
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1. Introduction 

Heavy fermion systems exhibit very interesting thermodynamic properties at low 
temperature. The Sommerfeld ratio of the specific heat to temperature at low T 
is very large in these materials indicating an effective mass for the electrons near 
the Fermi surface that can be orders of magnitude larger than the free electron 
mass, hence the term heavy fermions. In low-temperature experiments on certain 
heavy fermion alloys at critical doping the Sommerfeld ratio does not settle down 
to a constant, as it would for a conventional Fermi liquid, but continues to rise as 
the temperature is lowered [1]. This behavior has also been observed in systems at 
near-critical doping by tuning external parameters such as pressure or magnetic field 
[1, 2]. The non-Fermi-liquid behavior of the specific heat is attributed to strongly 
correlated physics, governed by an underlying quantum critical point (see [3, 4, 5] for 
reviews of heavy fermion systems) but a detailed understanding of quantum critical 
metals remains a key problem in theoretical physics [6] . 

The last few years have seen increased interest in developing dual gravity models 
for strongly interacting systems in condensed matter physics (see [7, 8, 9, 10, 11] 
for recent reviews). While the gravitational approach has found a number of in- 
teresting applications it should be emphasized that its ultimate relevance to real 
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condensed matter systems remains speculative. The original AdS/CFT conjecture 
[12], for which there is by now strong evidence, relates maximally supersymmetric 
3+1 dimensional Yang- Mills gauge theory with a large number of colors to super- 
gravity on a ten dimensional AdS5xS5 spacetime. In condensed matter physics one 
is instead interested in systems without supersymmetry where the validity of a grav- 
itational dual description is a more open question. Furthermore, the large number of 
colors plays a key role in the supersymmetric gauge theory context in justifying the 
use of large-scale classical geometry on the gravitational side of the duality and it 
is at present unclear what the corresponding formal limit would be for the strongly 
coupled condensed matter systems that one wants to model. These are important 
issues to settle but in the present paper we pursue a more modest goal of identifying 
yet another physical effect that can be phenomenologically modeled by a relatively 
simple dual gravitational system. 

Gravity duals have one or more additional spatial dimensions compared to the 
field theory in question. Finite temperature effects are studied by having a black hole 
in the higher-dimensional spacetime and when the black hole is electrically charged 
the Hawking temperature and the black hole charge provide competing energy scales 
that can lead to interesting dynamics. We will consider quantum critical points in 
three spatial dimensions that are characterized by a scaling law. 



where the dynamical critical exponent z is in general different from 1. We arc par- 
ticularly interested in quantum phase transitions in heavy fermion alloys where the 
metal goes from being paramagnetic to being antiferromagnetic as the level of doping 
or some other external control parameter is varied [1, 2]. Theoretical work based on 
a lattice Kondo model suggests that the associated quantum critical point exhibits 
anisotropic scaling with a non-universal dynamical critical exponent that can be fit- 
ted by comparison to experimental data [13]. In particular, it was found that z ~ 2.7 
for the critically doped heavy fermion alloy CeCus.gAuo.i, whose anomalous specific 
heat was reported in [1]. 

In [14] it was conjectured that a strongly coupled system at a fixed point governed 
by anisotropic scaling of this form could be modeled by a gravity theory in a so called 
Lifshitz spacetime with the metric,^ 



X ^ Xx , 



(1.1) 




(1.2) 



which is invariant under the transformation 



t yt, X ^ Xx, r -)■ - . 



r 



(1.3) 



See [15] for early work on gravitational backgrounds with anisotropic scaling. 
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Here L is a characteristic length scale, which will determine the magnitude of the 
negative cosmological constant in the gravity model, and {t, x, r) are dimensionless 
coordinates in the higher dimensional spacetime. The finite temperature physics of 
such a system is then encoded into the geometry of a black hole, which is asymptotic 
to the Lifshitz spacetime (1.2). 

In this paper we continue our study of asymptotically Lifshitz black holes from 
[16, 17] placing an emphasis on their thermodynamic properties. We compute the 
specific heat of charged asymptotically Lifshitz black holes as a function of tem- 
perature, keeping the charge density in the dual field theory fixed. At high T the 
thermodynamic behavior of the dual system is governed by the symmetries of the 
quantum critical point leading to a temperature dependence of the specific heat of 
the form C ~ T^/^ This is a direct consequence of scaling and we find the same 
result in the gravity dual using black hole thermodynamics. 

At low T, we observe a crossover in the black hole thermodynamics indicating 
non-trivial collective effects in the dual field theory. For z — 1 our system reduces 
to the thermodynamics of the well-known AdS-RN black branes. In this case, the 
specific heat divided by T goes to a constant value in the limit of low temperature, 
which is the same behavior as seen in a conventional Fermi liquid. This result was 
obtained previously in [18], where it was interpreted as evidence for the presence of a 
Fermi surface in the dual system. It was noted in [18], however, from a holographic 
computation of transport properties that the low-energy effective theory near the 
Fermi surface cannot be that of weakly interacting fermions. Non-Fermi-liquid be- 
havior has also been observed in spectral functions of probe fermions coupled to the 
z^l system [19, 20]. 

For a non-trivial dynamical critical exponent z > 1 the black brane thermody- 
namics calculation is more involved and we have to rely more heavily on numerical 
computations. Interestingly, we find qualitatively different behavior in this case. The 
Sommerfeld ratio C /T now continues to rise as we go to lower temperatures, which is 
indeed the behavior observed in critically doped heavy fermion systems. The ability 
of gravity models with Lifshitz scaling to capture this non-Fermi-liquid aspect of real 
quantum critical metals is interesting and is the main result of the present paper. 

2. The holographic model 

The application that we have in mind involves a dual field theory in three spatial 
dimensions and accordingly the gravity dual is defined on a 4+1 dimensional space- 
time to accommodate an emergent extra dimension. Most of our results carry over 
to other dimensions in a straightforward fashion. In particular, similar conclusions 
can be drawn about planar systems with d — 2. 

We work with a 4 -|- 1 dimensional version of the holographic model of Kachru 
et al. [14], as formulated by Taylor [21], and with a Maxwell gauge field added, as 
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in [17]. The bulk classical action for the gravitational sector consists of two parts, 

S — <S'Emstein-Maxwell + "S'Lifshitz ■ (2-1) 

The first term is the usual action of 4 + 1 dimensional Einstein-Maxwell gravity with 
a negative cosmological constant, 

'S'Einstein-Maxwell = J — 2k— -F^j^F'^"^ . (2.2) 

This is followed by a term involving a massive vector field, 

Suf.m. = -J d'x,/^ {jV'' + ^A^') ' (2-3) 

whose sole purpose is to provide backgrounds with anisotropic scaling of the form 
(1.3). We refer to this auxiliary vector field, which only couples to gravity, as the 
Lifshitz vector field. The original gravity model of [14] was 3-|-l-dimensional and 
anisotropic scaling was obtained by including a pair of coupled two- and three-form 
field strengths. Upon integrating out the three-form field strength, the remaining 
two-form becomes a field strength of a massive vector and in this form the model is 
easily extended to general dimensions [21]. 

The equations of motion obtained from the action (2.1) consist of the Einstein 
equations, the Maxwell equations, and field equations for the Lifshitz vector, 

Gi A „ rpMaxwell i rpLifshitz A\ 

ixv + ^^QfJiv — -I- +-'-nu : 

V^F'^^ = 0, (2.5) 

V^^J'"" = a'A" , (2.6) 
with the energy momentum tensors of the Maxwell and Lifshitz vector fields given 

by 

yMa^eii ^ ^^F,,F^^ - ^g,^F,„F'n, (2.7) 
^Lifshitz ^ 1 _ Ig.uJ^x.J^'n + ^{A,A. - lg,.AxA'). (2.8) 

The Lifshitz fixed point geometry (1.2) is a solution of the equations of motion 
provided the characteristic length scale L is related to the cosmological constant A 
via 

A = -^, (...) 



the mass of the Lifshitz vector field is fine-tuned to c = y/Sz/L, and the Lifshitz 
vector field has the background value 



A^W ^^^ ^h r', A,,^Ar^O. (2.10) 
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The Maxwell field vanishes in the Lifshitz background, = 0. 

Although we do not do it here, it is straightforward to couple matter to this 
system. Including a scalar field, for instance, leads to an instability at low T giving 
rise to holographic superconductors with Lifshitz scaling [17]. It is also interesting to 
include a coupling to a bulk Dirac spinor and obtain fermion spectral functions along 
the lines of [19, 20, 22]. We find that the spectral functions in the z > 1 theory exhibit 
many of the same features as in the asymptotically AdS bulk spacetime considered 
in these earlier works, including peaks at certain values of momentum and frequency, 
but there are also important differences. In particular, the would be quasiparticle 
peaks are less sharply defined than their z — 1 counterparts.^ 



3. Charged black branes 

In order to study finite temperature effects in the dual strongly coupled field theory 
we look for static black brane solutions of the equations of motion (2.4) - (2.6) 
which are asymptotic to the Lifshitz fixed point solution (1.2). We consider a 4+1 
dimensional metric of the form 

ds'' = L' [-r^'firfde + r'dx^ + ^-^dr^^ , (3.1) 
for which the non- vanishing components of the vielbein can be taken to be 

e'i^Lr^fir), ^ . . . ^ ^ Lr , = L ^ . (3.2) 

An asymptotically Lifshitz black brane with a non-degenerate horizon has a simple 
zero of both /(r)^ and 5'(r)~^ at the horizon, which we take to be at r = ro, and 
/(r), g{r) — >■ 1 as r — >■ oo. It is straightforward to generalize this ansatz to include 
black holes with a spherical horizon or topological black holes with a hyperbolic hori- 
zon but it is the fiat horizon case (3.1) that is of direct interest for the gravitational 
dual description of a strongly coupled 3 + 1 dimensional field theoretic system. 

In a static electrically charged black brane background the Maxwell gauge field 
and the Lifshitz vector can be chosen to be of the form 



AM^{a{r),0,...,0}; Am = \J '^^^j^{a{r),0, . . . ,0} , (3.3) 

where the corresponding coordinate frame components are given by ^4^ = c^Am and 
A, = e^Au. 

^Work in progress. 
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3.1 Field equations for static configurations 

The Maxwell equations and the equations of motion for the Lifshitz vector can be 
written in first-order form, 

r d 
J dr 
d 



r d 

dr ^ ' 



if a) = -za + g/3, 
0, 

—z a + z gb, 



(3.4) 
(3.5) 
(3.6) 



3r ga, 



(3.7) 



where /3 and b are defined via (3.4) and (3.6) respectively. The Maxwell equation 
(3.5) trivially integrates to 

/3(r) = ^, (3.8) 

where the integration constant p is proportional to the world-volume electric charge 
density of the black brane, which in turn is identified with the charge density in the 
dual field theory by the standard AdS/CFT dictionary [23, 24], as explained in detail 
in [25]. 

The Einstein equations reduce to the following pair of first-order differential 
equations 



rdg 
gdr 

f dr 



3_ 
3 

3 



iz-l 



2 V ; 2 4r6 



P 

^^6 



+ 2, 



-z-l. 



(3.9) 
(3.10) 



2 ^ ' 2 

The Lifshitz fixed point solution is given by a = p = and f = g = a = b = l. 

To bring out the underlying scale invariance of the model we rewrite the field 
equations using u — log(r/ro) instead of r and introduce the scale invariant ratio 
p = p/r^. The remaining Maxwell equation (3.4) and the equations of motion of the 
Lifshitz field (3.6) and (3.7) reduce to 



a + —a — —za -\- pe ^g . 
f 

d + —a — —za -\- zgb , 
h — —36 -|- "iga , 



(3.11) 

(3.12) 
(3.13) 



where ' = ^. The Einstein equations (3.9) and (3.10) can similarly be rewritten 



du ' 

9_ 
9 
I 
f 



9_ 
6 

6 



[z - 1) (3a' + zb'') + ^e"^" - {z" + 2^ + 9) 
{z - 1) (3a2 - zb"") - ^e"^" + {z^ + 2;^ + 9) 



+ 2, 
-z-\. 



(3.14) 
(3.15) 
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Using these variables, the field equations are manifestly invariant under the rescaling 
(1.3). For given z > 1, there is a one parameter family of black brane solutions, 
labelled by p. A neutral black brane has p = while the extremal limit is given by 

3.2 Exact solutions 

When z — 1 the terms involving a and b on the right hand side of the Einstein 
equations drop out. The equations then reduce to those of Einstein-Maxwell theory 
with a negative cosmological constant and one finds the AdS-RN solution describing 
charged black branes in 4 + 1 dimensional asymptotically AdS spacetime. 



f = l = (l-e-) (^l + e--ge-«) 



(3.16) 



A = ^ (1 - e-^") . (3.17) 

At z > 1 there is a corresponding family of charged asymptotically Lifshitz black 
branes. In general these solutions can only be obtained numerically, but it turns out 
that in the special case z — 6 an isolated exact solution can be found [17, 26], 

6=1, /2 = l = a' = l-e-^", At = ±V^Lr^(e3"-l). (3.18) 

We will primarily rely on numerical solutions in the following but the AdS-RN solu- 
tion and the exact z = 6 solution are useful for checking the numerics. 

3.3 Numerical solutions 

The family of black brane solutions at generic z > 1 can be mapped out using 
numerical techniques similar to those of [16]. The numerical integration is started 
near the black hole, with suitable boundary conditions for a regular non-degenerate 
horizon at -u = 0, and then proceeds outwards towards the asymptotic region. First 
of all, there should be a simple zero of and g^"^ at the horizon. Furthermore, the 
product fa must have a simple zero at the horizon in order for the gauge connection 
to be regular there. Similarly, the combination ga on the right hand side of (3.13) 
should be regular at the horizon and it then follows from (3.12) that b takes a finite 
value there. Putting all of this together we find that the near-horizon behavior can 
be parametrized as follows 

f{u) = V^(/o + M + + ...), (3.19) 
g{u) = ^{go + giu + g2U^ + ...), (3.20) 

a{u) — \/u{ao + aiu + a2U^ + ■■■), (3-21) 
a{u) — \/u{aQ -\- aiu + a2U^ + ■■■), (3.22) 
b{u) = 6o + hu + b2U^ + .... (3.23) 
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Inserting these near-horizon expansions into the equations of motion we obtain a 
two-parameter family of initial value data for any given 2;, parametrized for instance 
by p and 60 • It turns out, however, that these two parameters cannot be varied 
independently but are restricted by the condition that the metric and Lifshitz vector 
approach the Lifshitz fixed point solution (1.2) and (2.10) sufficiently rapidly as 
M — > 00 [16, 27, 28]. This means that for a given value of z the parameter p uniquely 
determines a charged black brane solution, up to overall scale. 

3.4 Asymptotic behavior at large u 

We now consider the asymptotic behavior of our fields as m — )■ 00. The first observa- 
tion is that the system of equations (3.11) - (3.15) only contains the combination j 
and therefore a uniform rescaling of / is a symmetry of the equations. The symmetry 
is fixed, however, in an asymptotically Lifshitz solution for which / — j-lasw— j-oo.^ 

We next observe that (3.15) can be used to eliminate / from (3.12) - (3.14) 
leaving a closed system of first-order equations for a, 6, and g. The large u behavior 
can be obtained by linearizing around the Lifshitz fixed point a = b = g = 1. The 
corresponding problem for d = 2 was discussed in [16, 27] and we omit the details 
here. At large u, solutions of the full non-linear system approach a linear combination 
of the eigenmodes of the linearized system plus a universal mode coming from source 
terms involving p^. The leading large u behavior of the eigenmodes is 0{e^'^) with 
the eigenvalues Aj G {—3 — z,\ (—3 — 2; ± V9-2^ — 26^; + 33) } while the universal 
mode falls off as 0(e~^") independent of z. 

Figure 1 shows the three eigenvalues as a A 
function of z. The eigenmode that belongs to , — — \ ^ 



finite energy. This mode is eliminated in our 

black hole solutions by fine-tuning the value Figure 1: Eigenvalues of the lin- 
of h at the black hole horizon. For z < 3 the earized problem for {5, 6, a} as a func- 
universal e~^" mode is sub-leading compared ^"^^ ° ^ or a — cs. 
to the two remaining eigenmodes of the lin- 
earized system but at 2; > 3 it dominates the asymptotic behavior. In the dividing 
2; = 3 case the linearized system is degenerate and the asymptotic behavior includes 
powers of u on top of the e~^" falloff. 

Finally, we turn to equation (3.11) for the Maxwell field. In the Lifshitz back- 

^When solving the equations numerically, it is convenient to initially put /g = 1 in the near- 
horizon expansion (3.19) when setting up the numerical integration and then rescale / at the end 
of the day so that / — 1 as u — 00. 



the largest eigenvalue is problematic [16, 27, 
28]. For 2; > 3 it is non- negative and a solu- 
tion that contains it fails to be asymptotically 
Lifshitz. For 1 < 2; < 3 it is a negative mode 
but the falloff at large u is too slow to give 




-4 
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ground with f — g — 1 it has the solution 



/X e + p M e 11 z — 3. 



The non- vanishing component of the vector potential in a coordinate frame is then 

(^/x + Lplog(^) if ^ = 3, 

with II — LjlrQ. Once again the qualitative behavior of solutions changes at 2; = 3. 
For low z values, in the range 1 < 2; < 3, the term involving the chemical potential 

/i dominates compared to the term involving the charge density p at large r, while 
at 2; > 3 it is the charge density term that is leading [29]. 

In general, we are not working with Lifshitz background itself but with solutions 
of the full non-linear system of equations that are only asymptotically Lifshitz. It 
turns out, however, the leading large u behavior of a carries over from the Lifshitz 
background to more general case as long as the value of z isn't too high.'' 



4. Black brane thermodynamics 

So far, we have set up the holographic model and considered charged black brane 
geometries that are conjectured to provide a dual description of a strongly coupled 
3+1 dimensional system near a quantum critical point with dynamical critical expo- 
nent z > 1. These black branes turn out to have interesting thermodynamics. At 
high temperature the thermodynamic behavior is governed by the underlying Lif- 
shitz symmetry, but collective effects come into play as the temperature is lowered 
and modify the thermodynamics. To see this, we use a combination of analytic and 
numerical arguments to compute the specific heat of the system as a function of tem- 
perature. Our starting point is the following expression for the Hawking temperature 
of a black brane, ^ 

Th = ^^-^, (4.1) 
47r go 

which is obtained in the standard way by requiring the Euclidean metric of the black 
brane to be smooth at the horizon at r = tq. The coefficients /o and go are taken 
from the near-horizon expansions (3.19) and (3.20) in the scale invariant variable u 



^At z>9 we expect non-linear effects to give rise to additional terms in (3.24) with a falloff in 
between that of the charge density and chemical potential terms. We will not keep track of such 
terms here since all systems of physical interest that we are aware of have z < 9. 

^Prom here on we fix the characteristic length scale as L = 1. Explicit factors of L can be 
reintroduced into the formulas by dimensional analysis. 
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and only depend on the value of p. The temperature can therefore be expressed as 
follows, 



H 



(4.2) 



with Fz{p) = ^ for different values of z shown in Figure 2. 

The specific heat at fixed volume in 
the boundary theory is then 



C = T 



dS 



T 



(dS/dro) 



(4.3) 




Figure 2: Temperature function Fz{p) for 
several z values. The z = 1 curve plots the 
exact result (4.10) while the z > 1 curves are 
obtained from numerical black hole solutions. 



dT (dT/dro) ' 

where 5* = is the Bekenstein-Hawking 
entropy density of the black brane and T 
in the boundary system is identified with 
the Hawking temperature (4.2). We have 
to make a choice whether to work at 
fixed charge density p or fixed chemi- 
cal potential p in the boundary theory 
when calculating thermodynamic quan- 
tities. For the field theory problem that 

we wish to model it is natural to keep the charge density fixed since the net density 
of charge carriers is given by the density of dopants, which is fixed in a given sample 
[29]. In this case 

d d 3p d 
drQ Otq tq dp 
and we find 

1 = . (4.5) 

If we instead were to keep the chemical potential p fixed then 



(4.4) 



d 

dro 



_d_ 

dro 



zp 
ro 



dp 
dp 



d_ 

dp 



(4.6) 



with dp/ dp obtained from the asymptotic behavior of the Maxwell field (3.24) in the 
numerical black brane solution, and then the remaining steps in the calculation of 
the specific heat are parallel to those at fixed charge density. 

The variables tq and p on the right hand side of equation (4.5) refer to the 4+1 
dimensional black brane geometry while it is T and p that have direct interpretation 
in the dual field theory. The latter variables are expressed in terms of the former 
via the definition p = p/r^ and equation (4.2) for the Hawking temperature. We can 
numerically invert these relations in order to express the specific heat (4.3) in terms 
T and p. Figure 3 shows the map between brane variables and physical variables for 
z = 2 and similar maps are obtained for other z values. 
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4.1 Scaling at high temperature 

Below, we present numerical results for 
C /T obtained via the above procedure for 
several values of z but some information 
can be obtained by analytic arguments. 
In particular, we can extract a scaling re- 
lation for the specific heat at high temper- 
ature. The argument is the same for any 
number of spatial dimensions d. Writing 
p — p/tq, the expression (4.1) for the tem- 
perature remains the same except that 
the detailed shape of the function Fz{p) 

depends on d. The limit of high temperature at fixed p corresponds to large Tq and 
p — )■ 0. The field equations depend on p^ in a smooth way so we expect Fz{p) to be a 
smooth even function of p (this is also evident from the graphs in Figure 2) and the 
denominator in the d dimensional version of equation (4.5) reduces to ^^^^(O). As a 
result the temperature and the entropy density, and their first derivatives, depend 
only on the overall scale ro in the high-temperature limit and the electric charge 
carried by the black brane does not affect the dynamics to leading order. 



Figure 3: Contour map expressing ro and 
p in terms of T and p sX z = 2. 



T 



An 



> 



(4.7) 



The high-temperature behavior of the specific heat then immediately follows from 
the general expression (4.3), 

C ~ ~ T'^A, (4.8) 

which reduces, in particular, to C ~ T^/^ when d = 3. It is straightforward to go a 
step further and integrate both sides of (4.3) with respect to T to obtain the following 
relation between energy and entropy of the system in the high-temperature limit. 



E 



d 



d + z 



(4.9) 



recovering the result found previously for electrically neutral black branes in [30] . The 
scaling behavior can be traced to the underlying Lifshitz symmetry of the quantum 
critical theory. It is easy to see that a generic statistical mechanical system in three 
spatial dimensions with a dispersion relation of the form co exhibits the same 
scaling behavior 



4.2 Low-temperature behavior 

At low temperature, the black brane thermodynamics exhibits interesting behavior 
due to collective effects in the dual field theory and the simple scahng that is seen 
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at high temperature no longer apphes. The behavior is quahtatively different in 
conformal systems with z = 1 as compared to Lifshitz systems with z > 1 and we 
consider these cases in turn. 

At 2; = 1 we have the exact AdS-RN black brane solution (3.17) for which the 
Hawking temperature is easily determined. One finds 

Fi(p) = 4-^, (4.10) 
and the specific heat at fixed p is given by 

^ = (4.11) 

At low temperature the charge on the black brane approaches the extremal limit, 
p — > -\/24 and the specific heat depends linearly on temperature, 

A weakly-coupled Fermi liquid has a specific heat that is linear in T at low temper- 
atures and it is interesting to see the same behavior emerge in the low-temperature 
limit of a black brane in Einstein-Maxwell theory without having made any reference 
to specific matter fields in the calculation. This result was obtained previously in [18] 
where it was taken as evidence for the existence of a Fermi surface in the dual field 
theory. Spectral functions of probe fermions coupled to the z = 1 system computed 
in [19, 20, 22] also strongly suggest the presence of a Fermi surface. Having a specific 
heat that is linear in T and a Fermi surface does not imply that the system consists of 
weakly coupled fermions. Indeed, the scaling behavior of excitations near the Fermi 
surface differs from that of a Landau Fermi liquid in the probe fermion analysis of 
[19, 20] and the zero-temperature conductivity computed in [18] was found to be 
diffusive rather than ballistic, suggesting that disorder plays a role. 

At z > 1 wc do not have explicit analytic black brane solutions except the 
isolated z = 6 solution (3.18) and our results for the specific heat in z > 1 systems 
are therefore based on numerical black brane solutions with non- vanishing Lifshitz 
vector field. The key quahtative difference compared to the z — 1 case can be seen 
in Figure 2. The function Fi{p) goes to zero in a hnear fashion in the extremal 
limit p — > -\/24 but for z > 1 both Fz{p) and its first derivative go to zero in the 
extremal limit p — )■ \/2(z'^ + 2z + 9). It then follows from equation (4.5) that the 
ratio C /T diverges in the T — )■ limit for z > 1. The growth in C/T towards lower 
temperatures is evident in the z > 1 curves in Figure 4.® This trend is in qualitative 
agreement with the measured specific heat of certain heavy fermion alloys [4]. The 
high T behavior of the specific heat curves in the figure, on the other hand, matches 
the scahng law (4.8). 

^Finite numerical precision limits how far the curves for 2; > 1 in Figure 4 can be extended 
towards low T. 
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Figure 4: Black brane specific heat divided by temperature for several z values calculated 
at fixed charged density p = 1- The z = \ curve is obtained from the exact result (4.11). 

5. Discussion 

In this paper we have employed a relatively simple holographic description of a 
strongly coupled 3+1 dimensional quantum critical point with asymmetric scaling to 
model the anomalous specific heat found at low temperature in many heavy fermion 
compounds. The calculation is performed entirely within the context of black hole 
thermodynamics and the nature of the low-lying spectrum of excitations is thus 
investigated without introducing any specific matter probes. It would be very in- 
teresting to complement the results obtained here by a study of spectral functions 
for probe fermions, generalizing the z = 1 work of [19, 20, 22], and a calculation 
of conductivities. This requires us to extend existing gravitational techniques for 
obtaining spectral functions and transport coefficients to models exhibiting Lifshitz 
scaling with non-trivial dynamical critical exponent z > 1, and our work in this 
direction is in progress. 

One limitation of the present work concerns the non-vanishing black hole area 
in the extremal limit, which corresponds to having non- vanishing entropy at zero 
temperature in the dual system. This can for instance be remedied by coupling the 
system to a scalar field and include the back-reaction to the scalar hair carried by 
the black hole at very low temperatures. In this case the area of the black hole 
horizon shrinks as the temperature is lowered and the system no longer has finite 
zero-temperature entropy. This was shown in [31] for asymptotically AdS black holes 
with hair and we have obtained analogous results in our z > 1 model with scalar 
matter included [32]. The scalar hair corresponds to having a charged superfluid 
condensate in the dual theory and the instability to developing superfluidity bears out 
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the general expectation that scale invariant quantum critical matter is never found 
as the true ground state of a system. A further study of the interplay between the 
heavy fermion physics and superfluidity in these holographic models is an interesting 
avenue for further work. 

We find it intriguing that experimentally measured deviations from Fermi liquid 
behavior in critically doped heavy fermion alloys can be qualitatively reproduced by 
holographic duals with Lifshitz symmetry. The curves for z > 1 in Figure 4 suggest 
a power law for C /T at low temperature rather than a logarithm. The experimental 
data vary from one material to another and both logarithmic and power law fits are 
used [4]. The fact that our curves for z > 1 approach straight lines at low T should 
not be over interpreted but it is encouraging to see the right trend coming from 
simple gravitational models with Lifshitz scaling. 
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